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(54) Error coding method 



(57) A decimated and interleaved multiplication ficients determine the multiplication table content and 
table for finite fields as is useful in Reed-Solomon ordering, 
encoding computations. The generator polynomial coef- 
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' S^r^Lr JTrJl^^- '!lf *^ T^^^^ *° ^ «^ d««*«=es. and. more particularly, to 

error correcbon coding in such electronic devices. k « 

BACKGROUND TO THE INVENTION 

^^.r^^^n^^I^'^'^^°''■^ ^^^"^ correction coding in o«ler to o/ercome 

errors ansmg from the transmission oi storage medium. Fon«ard error-correction coding (FEC) systems add redun- 
T^rlf "Ik «> *at ii^e receiver can detect and correct errors using on y the rSiS S^life 

eliminates the need Ibr the receiver to send requests tor retransmission to the transmitter 

SS?-^ One of the more popular eror correction code types is Reed-Solomon code. Reed-Solomon codes areUock 
codes w.^ maximumdistance separation and are highly efficient in their use of redundancy. The most appealing aspect 
SJlTSi^Tr "^nJ! °f ^^^^ decoding algorithms. See for example. Wicke^uid Bha^ava 

EdsO Reed-Solomon Codes and Their Applications (IEEE Press. Piscataway, New Jersey. 1994) 
iSS^Li^Sli'V level View Of a Reed-Solomon coding system. Ablockof information bits. I. is encoded 

into a codewoKj C w^hich is a larger block that contains the original information and redundant bits. After transmission 

erat^ I from C+E and r v^ll equal I if the number of errors bits in E is within the correction capabilities of the coda 
S??rth*,S!^l * Ti^^*^ description of Reed-Solomon coding. In particular, a Wod. of bk information 
^^^^S^^^t ^'^ *^ represented by a symbol, producing a btock of k infor- 

mahon ^mitx>ls for coding. The encoder operates on the block of k infomiation symbols to produce a block of n code- 
!I!^«f ^ ^ ttie original information in some form as well as redundancy The code can be designed so 

that the redundancy is used for error detection only, error correction only, or a combinatton of some error detec^on and 

channel. The receiver fiont-end produces a btock of bn bits that might be corrupted, depending upon the amount of 
^2? ^ t"" ^ "^''"^ ^ "<=«k °* " «y"*<"s and preceded with'tSe decSrSlg 

SJ,^ tTl^^f "^^ ^"'^^'^ ^y"*"'^' « ha^d<lecision decoder can reliably recover th! 

Tl^, J'T*'"" ^"^ '"P"* ""^ ^ P"°« tor the added redundancy is the inaease in the number 

i2??Li«ir«?H ^'^'S'" r*^"^ essentially maps k information symbols with the symbols as elements of a finite 
leM (Galo s leld or OF) wrth a power of 2 number of elements into n symbols which are GF elements from the sane 

mTr^FSwrra?^^^^^^ 

[MOT] Noreystematic Reed-Solomon encoding produces codewords by distributing the information and redun- 
dancy across the entire codeword of n symbols according to the coding algorithm. Systematic Reed-Solomon encod- 
ing, on ttie other hand, forms codewords by concatenating the k information symbols ith n-k pari^ sySoTv^S^ 
S^SS nS'"" '° ~f "9f 3°rithms. The additional n-k parity symbols contain the rSiSaSrn^IiWthS 
1^^. ^ T"'* K '"^y transmitted k information symbols. In particular, with receiver soft ded- 

Zt^u^to^MT^^'T *° "^"^ ^ symbols and detect s erased symbols provided aen-s is at 

m^ equal to n^.. Note that values such as n = 204 and k = 1 88 with the GF being GF(28) (the finite f iekl with 256 ele- 

r2M I8?r«2 « ^ode for high speed modems and would be called a 

(204, 1 88) code. This code can correct 8 error symbols per 204-symbal codeword 

S^!L,H * advantageous because the information component of the received 

JSI^ ^ infomiation. The last n-k symbols must be computed from the information symbols 
IrfrSJL™^^ "f*^ '^'^ be conputed from the information symbols using methods based on the arithmetic 

t^J^sZ coeff'oents are GF elements with the elements representing groups of bits. The information, par- 

S^n^^fT n-XTfA"**^ ^ polynomials l(x). P{x). and C(x). respectively For systematic Reed-Solomon 
encoding C(x) = x l(x) + P(x) with P(x) the remainder from the polynomial division of x"-'*l(x) by G(x) Gfx) is the 

ITxr->.r"x-VG ^-"'^^"V'^pf -""^ Polynomi-^ 'if dtS'^n'J 
«w-x +«n-k- ix +G^k-2X +GiX + Go,soP(x)hasdegreeatmostn-k-1.i(x)isaDolvnomial6f 

1 With coefficients being the n codeword symbols. ^ ^ ^« ..w 
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[001 0] Most popular architectures that implement polynomial division for systematic Reed-Solomon encoding com- 
prise feedback shift registers that are composed of delay elements, GF element nnultipliers, and GF element adders as 
shown in Figure 3 The delay elements D are initialized with zero symbol values. The information symbols are shifted 
into the register one at a time, highest order element first (l^.i). During each dock cycle of the register, the GF element 
5 held in the last delay element (leftmost] is fed back to n-k multipliers that compute the product of the feedback element 
with the feedback register multiplier elements through Gn.|c.v Because the finite field has a power of 2 number of 
elements, subtraction and addition are the same operation. 

[001 1 ] At each stage of the feedback register the products are added to the stored elements in the previous stage 
and the result is stored in the following stage. After clocking the register n times the elements stored in the delay ele- 
10 ments D are the remainder of the division, or the parity elements that constitute the coefficients of the parity polynomial 
P(x). Figure 4 shows a simplified feedback shift register that uses a pre-shifted l(x) to compute the remainder in only n- 
k clock cycles of the register. 

[001 2] It is important to understand that he architectures shown in Figures 3-4 evolved from a desire to efficiently 
implement Reed-Solomon encoders with circuit elements, in a typical encoder design, the three types of circuit ete- 

15 ments (delays. GF adders, and GF multipliers) are individually optimized and then put together to perform the desired 
remainder computation operation. This type of encoder architecture can be emulated on a general purpose digital sig- 
nal processing (DSP) platform. However, while either the GF multiply or the GF add can be implemented efficiently 
(depending upon the particular GF representation used), they cannot both be Implemented efficiently simultaneously. 
For example, one particular representation allows GF adds to be computed with a simple exclusive-OR operation of the 

20 binary components of the two elements. In general, this can be implemer^ed in one cyde of a DSP. However, for this 
same GF element representation, the GF multiply requires a large number of cycles to compute. 
[001 3] A GF multiplication table can be employed to reduce the number of cydes required to multiply two GF ele- 
ments. However, a GF multiplication table can require a large amount of memory, and memory lookup to determine the 
product of two GF elements can also be somewhat time consuming. 

25 

SUMMARY OF THE INVENTION 

[001 4] The present application cfisdoses a simplified finite field (Galois field or GF) multiplication table which is dec- 
imated and interieaved. 

30 [001 5] The present application further disdoses a method of performing polynomial division with f irnte field coeffi- 
cients. The method provides a finite field multiplication table with entries consisting of multiples of the coefficients of a 
divisor polynomial and ordered according to the coeff idents of said divisor polynomial. Partial quotients and reminders 
are then iterativeiy computed. 

[001 6] This has the advantages of allowing general purpose digital signal processors (DSPs) to efficiently perform 
35 Reed-Solomon encoding and thereby eliminates the need for specialized feedback shift register circuitry. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0017] Specific embodiments of the present invention will now be described in further details, by way of example, 
40 with reference to the accompanying drawings in which; 

Figures 1-2 schematically illustrate Reed-Solomon coding. 
Figures 3-4 show shift registres for polynomial division. 

Figure 5 illustrates a GF(256) multiplication tattle; and 
45 Figure 6 shows a preferred embodiment GF(256) multiplication table. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0018] Figures 5-6 illustrate the preferred embodiment GF(2^) element multiplication lookup table construction by 
so selection (dedmation) of columns from a full multiplication lookup table and reordering (interleaving) for efficient 
access. The coefficients of the Reed-Solomon code generator polynomial as elements of the GF indicate the columns 
to be retained, and the memory locations used for the emulation of the delays sets the ordering of the columns. Note 
for a (204,188) code, the code generator polynomial G(x) has degree 16i the information polynomial l(x) has degree 
187. and the parity polynomial P(x) has degree 15; so the polynomial division has 188 steps with each step requiring 
55 16 multiplications and 16 additions. GF coeff ident polynomial division 

[0019] The preferred ennbodiments emulate the feedt)ack shift register of Figure 4 with a general purpose proces- 
sor; thus first consider the shift register operation. The polynomial division of x^*^l(x) by G(x) to yield P(x) proceeds in 
steps with each step adding the next lower power of x term to the quotient and leaving a remainder of one lower power 
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Of X. The last step generates the constant term of the quotient and yields the coefficients of thermal remainder P(x) in 
the delays. Each step corresponds to one clock cyde in the feedback register of Rgure 4. Now: 



Q(x) = x-^^ Gn.Mx"-''-^ + ... + Go and l(x) = I ,.,x'^^ + 1 ...x'^'* + I ..gX ■'"^ + ... + 

Thus: 



x"-'i(x) » \^.,x'^' + i,.2x"-'' + 1 ..gx"-^' + ... + lox" ■^ 

10 And the first step of the division gives a first quotient term l^-ix^"^ and a fiist remainder equal to: 
("k-2- G„-k-ilM)x""' + Ok^-Qn.k-2" k-Ox"^ + ... + (Imk„.k)-GoI k-Ox**'' "" 

+ "k.2-(n.R)X'-^ + lk-3^n.k)X'"^ + ... + lox""'. 



During the first clock cycle the feedback shift register of Figure 4 shifts in I,,.,, computes the products Qlk i for j = 0 to 

n-k-1 . and^es them in the delays 0. 1 n-k-1 . The first dock shifted in 1^., but the terms Ik.2. Ik-a. ... . lo. have not 

yet been shifted in. The second step then gives a second quotient term: 



'k-2-Givk-i'k-i)*''^ 



And a second remainder 



('k-3- Gn-k-2'k-1-("k-2- <3n-k-1«k-l)Qn-k-l)x""^ + ('k-*' G„.K.3lK-r(>k-2- G„.k.ilfc.l) 
Gn.k-2)x'^ + ... + (l|^Hn^O-Go"k-i - (lk-2- GMt-ilk-OGOx"-' + (l,^.(^^ - 
('k-2- G„.,^ll K.,)Go)x''-'^ + Ik-S-Cn-lOX*"' + ... + loX"*. 



Dunng the second clock cycle the shift register shifts in 1^.2. computes the terms (l^.g- Gn-k 1 1k i)Gi for j = 0 n-k-1 
subiracls (adds) them (except for j = 0) from the terms Gj., I^.i . and stores in the delays. Note that the terms (l',- G^ J 
ilk.i)Gj, are computed by first subtracting (adding) the product Gn-k-i ik-i (stored in delay n-k-1 by the previous clock 
cyde) from the shifted in Ik-a and then multiplying the result by G in each of the multipliere. Also, the Q . L i term to add 
to this product were stored in the acyacent delay from the first ckx* cyde and is shifted in for the adcGtion' and storage 
Again, the terms Ik-a lo have not yet been shifted in. 

PMKO] Sinralarly. the successive division steps buikl up the remainders until the kth and last step has the remainder 
P(x) coefficients Pn^^., Pq in the delays n-k-1 0. 

[0021] The preferred embodiments rely on the following analysis of the foregoing polynomial division On the ith 
step the symbol. M„^.^ . stored In the n-k-1 memory Is subtracted (added) from the syirtbol lk.| shifted in. and the result 

IS multiplied by eadi of the generator polynomial coefficients, Qn-k-i Gq. Thus if the symbol lk.| - M„ k i is denoted 

a where a is the primitive element of the QF. then the n-k multiplications are a" k i. Ou.^ 9 and o"" Gn 
Hence, the multiplicands are all the same, namely, a"". 

[0022] Now, consider implementing these multiplcations by a lookup table. The element a" is used to index one of 
the rows of the multiplication table and the element G, indexes one of the columns of the table. The consecutive multi- 
plications executed during the do* cycle will require the same row to be accessed each time, the columns indexed cor- 
respond to elements Gn-k-i Gq. Therefore, the products produced by the consecutive multiplies are extracted from 

a given row and vanous columns of the table. In general, the column indices will not be consecutively oidered For any 
clock cycle of the register the row index can change, but the column indexing pattern will remain the same. 
[0023] The amount of memory required to store the multiplk»tion table can be reduced by eliminating the columns 
that are never indexed. The remaining columns in the multipication table correspond to the elements Gp k 1 Gq 
Thus the memory size to store the new multiplication table with only these columns is 2" by n-k where GF is the finite 
field GF(2'^. For example, with GF(28). the fuH table would be 256 by 256 elements, but with the (204.188) code the 
55 table would be only 256 by 16 elements. 

[0024] The number of DSP cydes required to implement the GF multiplies by table lookup can be reduced by oider- 

ing the columns so that they conrespond to the consecutive elements Gn.k.i Go (or Go. G, Gn.k.i. depending on 

the ordering of the memory locations emulating the delays). After ordering columns in this manner, multiplication by a"" 
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is implemented by accessing the appropriate row and then reading out the consecutive n-k elements of that row. 
[0025] The decimated multiplication table implementation can be further simplified by collapsing rt into a one- 
dimensional structure where the a"' element determines the offset into the one-dimensional memory array. Then, n-k 
consecutive elements are read from consecutive memory locations in the an^ay. In the case where each element is 
5 smaller than the accessible memory element size, several products can be accessed simultaneously. For example, if 
each product is an 8-bit GF element, and the DSP can access 32-bit words from memory, then four products can be 
read in one memory access. 

[0026] Figure 5 illustrates the conceptual manner in which the GF(2®) products of interest are extracted from a full- 
sized GF multiplication table and Figure 6 shows them placed into the decimated and interleaved GF multiplication 

10 table. In this example, each GF element is represented by 8 bits. Therefore the full multiplication table contains 256 by 
256 elements. The first row is trivially Os, and so are the first n-k elements of the new array. Starting with the second 
row (for element a° = 1). product element a°Go is placed in the new array, element a^G^ is placed in the second loca- 
tion, and so forth up to element a^Gn-k-i- Thus, the second n-k elements are the products of a° and the ordered ele- 
ments (coefficients) of the generator polynomial G(x). The products of a** and the ordered elements of G(x} are placed 

15 In the next n-k array locations, and so forth. During operatton, result of multiplying of ct} by the ordered elements of G(x) 
is obtained by accessing the n-k elements in the memory array starting at the index j(n-k). 

[0027] Figure 6 illustrates the decimated and interleaved GF multiplication table of 256(n-k) elements. An 8-bit input 

or multiplicand IsO or d for j=0. 1,2.. ...254. The shift register multipliers (the other multiplicands) are Gq. G-j Gp^.v 

The products are stored as an interleaved table so that all products for any input and the shift-register multipliers can 
20 be looked up in the appropriate order within one block. This makes it possible to perform multiple multiplications in one 
lookup cyde. For instance, a 32-bit load can produce the products for four 8-bit GF multiplications. This Is highly efficient 
for DSP implementation with wide memory busses. 

[0028] The preferred embodiments can be modified in various ways while retaining the feature of a multiplication 
table decimated to the Reed-Solomon code generator polynomial coefficient and interleaved for sequential access. 
25 [0029] For example, the finite field may have a different size; fhs sequential access may be adapted to another 
method of polynomial division; 

Claims 

30 1 . A method of performing polynomial division with finite field coefficients, which method comprising: 

providing a finite field multiplication table with entries consisting of multiples of the coefficients of a divisor pol- 
ynomial and ordered according to said coefficients of said divisor polynomial; and 
Iteratively computing partial quotients and remainders. 
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